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Abstract 

We determine the jet vertex for Mueller-Navelet jets and forward jets in the 
small-cone approximation for two particular choices of jet algoritms: the kt algo¬ 
rithm and the cone algorithm. These choices are motivated by the extensive use 
of such algorithms in the phenomenology of jets. The differences with the original 
calculations of the small-cone jet vertex by Ivanov and Papa, which is found to be 
equivalent to a formerly algorithm proposed by Furman, are shown at both analytic 
and numerical level, and turn out to be sizeable. A detailed numerical study of 
the error introduced by the small-cone approximation is also presented, for various 
observables of phenomenological interest. For values of the jet “radius” R = 0.5, 
the use of the small-cone approximation amounts to an error of about 5% at the 
level of cross section, while it reduces to less than 2% for ratios of distributions such 
as those involved in the measure of the azimuthal decorrelation of dijets. 


1 Introduction 

The identification of high-energy (Regge) dynamics in QCD processes has been long since 
theoretically investigated. With the advent of high-energy colliders like HERA, Tevatron 
and LHC such studies have become possible at experimental level too. The processes 
which are expected to be more sensitive to this peculiar dynamical regime, where the 
center-of-mass energy ^/s is much larger than all hard scales involved in the scattering, 
are the so-called Mueller-Navelet (MN) jets at hadron-hadron colliders [T] and forward 
jets at electron-hadron colliders [2]. Both of them are dehned by the presence of QCD 
jets at large rapidity, accompanied by any hadronic activity which is inclusively collected 
in the central region. 

Such processes can be theoretically described by factorization formulae which involve 
several ingredients: the partonic distribution functions (PDFs) of the incoming hadron(s), 
the gluon Green’s function (GGF) describing the high-energy dynamics of emitted and 
exchanged partons — mostly (reggeized) gluons — and hnally the so-called jet vertices, 
describing the production of a forward jet from the interaction of one incoming parton 
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and a reggeized gluon. In the case of an incoming electron, an additional quantity, the 
photon impact factor, has to be considered too. 

At present, all ingredients are known at next-to-leading level in the respective parame¬ 
ters: the PDFs which resum logarithms of collinear type, the GGF resumming logarithms 
of the energy, and the jet vertices (and impact-factors) which are computed at hnite 
perturbative order. 

Focusing for dehniteness on MN jets in hadron-hadron collision, the pro cess-dependent 
part of the cross section is represented by the jet vertex [3], which depends on the jet 
variables and the actual jet algorithm. By following the work of Ivanov and Papa i , in 
this paper we reconsider the computation of the jet vertex in the small-cone approximation 
(SGA), namely for jets whose extension in the rapidity-azimuthal angle (?/, 0)-plane is 
small. In particular we apply their method to determine the analytic expressions of the 
jet vertex for two particular choices of jet algorithms: the kt algorithri% [S] and the cone 
algorithm [B]. These are the mostly used algorithms in modern jet phenomenology, in 
particular the kt one. On the contrary, the algorithm used in |1] can be traced back to 
the one considered by Furman [7] in early studies of QGD radiation, but not used for 
practical purposes anymore, being infra-red unsafe. 

The aim of our work is twofold: on one hand we want to give a precise estimate of the 
error introduced by the small-cone approximation in the description of QGD observables at 
high energies, i.e., at large rapidities, so as to possibly justify its use in phenomenological 
analyses. On the other hand, we want to estimate the differences occurring by choosing 
different jet algorithms for the same process. The jet algorithm dependence has already 
been studied in the past, and a detailed analysis of the cone and kt algorithms in the 
SGA was presented in ref. [S] in the context of collinear factorization. Here we carry out 
a similar analysis in the framework of high-energy (kt-dependent) factorization, the basic 
tool for the description of the Regge regime in perturbative QGD. 

One should keep in mind that the small-cone expressions are fully analytic (before 
their convolution with the PDFs) and compact, and allow a simple implementation in 
numerical codes that run much faster than those with the exact jet vertices. For this 
reason the SGA jet vertices have already been used PDTO] in quantitative comparison 
with available data m- However, the experimental results were extracted by clustering 
jets with the kt algorithm, while the SGA jet vertices used for the theoretical calculation 
were those obtained with the Furman algorithm by Ivanov and Papa (FIP). Also the small- 
cone analysis performed in ref. [12] compared calculations with the exact jet vertices in 
the kt algorithm versus the small-cone ones in the FIP algorithm. 

Starting from these premises, and after reviewing in sec. |2]the theoretical setup for the 
description of Mueller-Navelet jets in terms of the collinear and high-energy factorization 
formulae, in sec. |3]we discuss in detail the differences among the three (kt, cone and FIP) 
jet algorithms in the relevant case of two near particles and we determine the kinematical 
conhgurations in the limit of small jet “radius” R. 

We then derive the small-cone jet vertices for the kt and cone algorithms in sec. S] by 
computing their differences with respect to the FIP algorithm induced by the different 
kinematical conditions. 

In sec. El we perform a numerical study in order to asses the quantitative difference 

^Here kt algorithm denotes the whole class of clustering algorithms based on ref. 0, which may 
differ in the details of the recombination scheme and of the resolution variable, like the anti-kt and the 
Cambrigde/Aachen versions. 
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between the exact and small-cone jet vertex in the kt algorithm, by comparing the vertices 
themselves as well as a typical differential cross section for MN jets and some angular 
coefficients measuring the azimuthal decorrelation between the jets. In addition we also 
determine the discrepancies of the same quantities induced by a different choice of the jet 
algorithm. 

We discuss the results in sec. El where we conclude that the wrong choice of algo¬ 
rithm causes sizeable errors on the predictions, while the SCA within the same algorithm 
provides a good approximation to the exact quantities and can therefore be used as a 
valuable tool for a quantitative description of MN and forward jets. 


2 Theoretical setup 

2.1 Factorization 

The process we are considering was suggested long ago by Mueller and Navelet pQ in 
order to study the high-energy behaviour of QCD. It is generated by the collision of two 
hadrons Ha,b — typically (anti)protons — and is characterized by the detection in the 
hnal state of two hard jets Ji ^2 with large rapidity separation; 

Ha + Hb^ J 1 + J 2 + X (1) 


where X represents any additional emission. Each jet Jj represents a cluster of particles 
grouped together according to some given jet algorithm and is described by 3 variables: 
the rapidity ?/j, the transverse energy Ei = |fcjj| and the azimuthal angle = arg(fcj_j), 
kj^i being the i-th jet transverse momentum. 

The kinematical region where one expects the high-energy QCD dynamics to play an 
important role is given by 


S={pa+ PbY > Ef ~ , \Y\ = \y^ - y^l 1 , (2) 

where the condition of hard jets {Ef S> Aq^^) is imposed for the applicability of pertur¬ 
bation theory. 

In the leading twist approximation, i.e., up to power suppressed correction in the 
hard scale parameter Aq^^/E^ 1, the hadronic cross section a can be factorized in 

the (longitudinal momentum fraction) convolution of two partonic distribution functions 
(PDFs) fa/nix) and a partonic cross section a 


daABjs) 

dJidJs 



da;ida;2 fa/A{xi)fb/B{x2) 


daab{xiX2s) 

dJi dJ 2 


dJi = dyi dEi d^* , (3) 


where a,b E {q, g} denote the parton flavours (quark or gluon) and Xi the partonic 
momentum fractions w.r.t. their parent hadrons. In turn, in the high-energy Regge regime 
we are considering, the partonic cross section for jet production can be factorized in a 
(transverse momentum) convolution of (process dependent) jet vertices V and a universal 
factor G called gluon Green’s function (GGF) 


ddab{XiX2s) 

d Ji d J 2 


j d^fei d^A;2 Va{xi, fep Ji)G{xiX2S, ki, k2)Vb{x2, ^2; J2) , 


(4) 
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Figure 1: Diagrammatic representation of the collinear and high-energy factorization for¬ 
mula for Mueller-Navelet jet production: fa^ represents the parton densities, Va^ the jet 
vertices and G the gluon Green’s function. 

where ki denotes the (reggeized) gluon transverse momentum flowing out from the GGF 
and entering the jet vertex, while s = X 1 X 2 S is the center-of-mass energy squared of the 
partonic subsystem. The overall factorization structure is depicted in fig. [H 

The PDFs are non-perturbative objects that depend also on the renormalization and 
factorization scales fiR, pip. While pip is introduced in the renormalization of UV diver¬ 
gencies, pip enters in the treatment of the IR collinear divergencies which are absorbed by 
the PDFs. The pip dependence of the PDFs is governed by the DGLAP equations [13], 
and their evolution kernels (splitting functions) are known at next-to-next-to-leading or¬ 
der (NNLO). By means of global fits they have been determined in a wide range of the 
{x, p.p) plane. 

The GGF is the central object in high-energy QGD, in that it resums the log(s) to all 
orders in perturbation theory. It obeys the BFKL equation m 

uiGJki, k,j) = S^(ki + fej) + j ifk K(ki, k)GJ,k, k,,) (5) 

/ ioo j / s \ 

— - GUk„k2), (6) 

-ioo27ri \soJ 

where is the Mellin transform of G(s), defined in terms of an arbitrary energy scale Sq. 
The BFKL kernel K = is known in next-to-leading-logarithmic (NLL) 

approximation. The coefficient depends on the choice of both fip and Sq. 

Finally, the jet vertices are perturbative finite objects without energy (\/i) depen¬ 
dence, and are known in NLO approximation: V = They depend on 

the jet variables y,E,<p and also on the arbitrary scales pip, pip, Sq, in such a way that 
the hadronic cross section be independent of those scales up to NLL terms, i.e., the scale 
dependence is present only in the terms of relative order a‘^{as log(s))"'. 

It is apparent that the determination of the Mueller-Navelet (MN) jet cross section 
involves quite a number of integrals, both in the factorization formulae fl3E|) . and in the 
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determination of G and V, as can be checked from their explicit expressions [3]. Further¬ 
more, when comparing the theoretical predictions with experiments, one needs integrated 
cross sections in some of the jet variables, in order to comply with the experimental 
binning. Such integrations are mostly done numerically, and this would require a large 
amount of computing resources or time for reaching a precision at the level of 1%. 

In order to cope with such a problem, two techniques can be exploited so as to reduce 
the computing time and improve convergence: 

• to project the GGF and jet vertices on a complete set of functions which respects 
the symmetries of the process (e.g., azimuthal invariance); 

• to use an approximated and simpler version of the jet vertices. 

These methods are often used in BFKL phenomenological analyses, and we shall illustrate 
them in the following subsections: the former in order to set up the theoretical framework 
and the main notations; the latter in order to introduce the main subject of this paper. 

2.2 Representation in Mellin space 

The hrst method to reduce computing time does not involve any approximation, at least 
at the NLL level of accuracy we are working with. It is better illustrated in the LL 
approximation, where the strong coupling ag is hxed and thus the BFKL kernel and jet 
vertices, in addition of being invariant under azimuthal rotations, are also scale invariant!^ 
In this case, a Fourier-Mellin transform diagonalizes the transverse integrations and the 
ensuing expressions are considerably simpler to evaluate. 

One can proceed in this way: hrst of all, let’s exploit the fact that, in the LL ap¬ 
proximation, each partonic momentum fraction coincides with the corresponding jet’s 
longitudinal momentum fraction (because of a (5(x — xj) in 1/*'°^): 

X = xj = (-I-, — for jet 1, 2) , (7) 

so that s = X 1 X 2 S = E 1 E 2 , where Y = \yi — 7 / 2 ! is the rapidity distance between the 

jets. If we adopt the convenient and natural choice of Sq = E 1 E 2 as energy scale, the 
Green’s function ([H]) is now independent of the partonic momentum fractions: 

/ ioo j 

— e-"'G^(fci, k 2 ) = G{Y, fci, k 2 ) , (so = E,E 2 ) (8) 

-ioo 

Secondly, we introduce the impact factor by integrating a jet vertex with the corre¬ 
sponding parton density: 

<hA(fc; = dx fa/A{x)Va{x, k] J) . (9) 

a -^0 

SO that we can rewrite the factorization formula in the form 

= j d^fci d^k2 <!>A{ki)G{Y, fci, k2)d>B{k2) . (10) 

^Strictly speaking, in our notations they are homogeneous functions of the transverse momenta fc, kj. 
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At this point we project kernel and vertices onto the eigenfunctions of the LL BFKL 
kernel 

Snuik) = , (n e Z, z/ G M) , (11) 




TT 


satisfying the completeness relation 

/ OO 

di/ £„,{k)£:,{k') = - k') 

■OO 


( 12 ) 


n£7j ' 


(0) 


and providing the LL eigenvalue function y- 


= / d^k' K^^\k, k')8Uk') = X^'nl^nuik) 


X^nl = 2^/^(l) - ^ 


1+77, 


+ iz7 -Ip 


1 + n 


— IV 


(13) 

(14) 


Finally, by inserting a completeness flT^ between each pair of factors in ffTOj) . we arrive 
at the convenient expression for the differential cross section 


da 


AB 


j dv ^Anu G^y)^ 


B nu 5 


(15) 


d J, d Jj 

where in the last equality we have used the Fourier-Mellin transforms 

j d^fc <F(fc; J)Sn.{k) = ^nu{J) (16) 

j d^k d^k' S:Ak)G{Y, k, fc')^nv(fc') = GUyK-^rS^n'S{v - v') . (17) 

The delta functions on the r.h.s. of eq. fll7p are just a consequence of the azimuthal- and 
scale-invariance of the kernel, and allow us to trade two bidimensional integrals for a sum 
and a simple integral. 

The azimuthal correlation of the MN jets is usually measured by means of the Fourier 
coefficients (m G Z) 


-271 


Gm{Ei,yi]E 2 ,y 2 ) = / d 0 id 02 cos (m( 0 i- 02 - 77 ))- 


da 


(18) 


d Ji d J 2 

Because of the azimuthal- and scale-invariance of the vertices, it is easy to show that 

E, 0) = , d’nu = (19) 

where T is dimensionless jl thus obtaining a factorization formula with just one integration: 


= 


277 

El E 2 


dv 


El 


2 \ 


d^Amu Gmu{Y) 


B mu 


( 20 ) 


Such a structure is preserved in the NLL approximation too, provided the impact 
factors and GGF are suitably modihed in order to take into account the loss of scale- 
invariance due to the renormalization procedure and to the factorization of collinear sin¬ 
gularities, which translates in a dependence on \og{E/fiji), log{E /and log(i7/yGo) of 


^We shall sometimes refer to the reduced impact factor as jet vertex. 
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the GGF and impact factors. The expression of the GGF in Mellin-space is very simply 
expressed in terms of the eigenvalue Xnv of the BFKL kernel; 

Gny{Y) = e^^"" , Xnu = d.Xnj + (21) 

However, the computation of ^nv at NLL level involves several integrations. The main 
advantage of this procedure is that such integrations can be done once and for all for each 
set of one-jet variables. Nevertheless, such computations can still be rather lengthy, and 
the use of an approximate expression of the impact factors — to be described in the next 
subsection — turns out to be very convenient. 

Let us conclude this section by noticing that the factorization formula (|20|) is very 
useful also in the case of cross-section integrated in jet energies. In fact, a double integral 
in El and E 2 factorizes into the product of simple integrals of the impact factors (provided 
the integration domain D = Ii x I 2 can be factorized into the cartesian product of two 
one-dimensional sets) 


Cm{yi,y2)= / dEidE2C„ 

JD=hxl2 

= {27rf j du GmAY) 

^ (27r)2 [ duGUY) ^AmAyi) ^BmAy2) 


( 22 ) 


[ dEi 


[ dE 2 {Ely’^-^dlBm. 

Jh 


Jh 


In this case, the integrated impact factors can be computed independently and stored in 
suitable grids, thus reducing a lot the computational effort of the phase-space integration. 
In the general case D A h x h, however, the expression fl20|) has to be numerically 
integrated in energy (and possibly in rapidity), and a suitable approximation (like the 
SGA) could be a very valuable tool to diminish the computing demand. 


2.3 Small-cone approximation 

In order to study the behaviour of the jet vertex for small values of the “radius” R — 
to be precisely dehned later on — and possibly to speed up the computation of the jet 
impact factor and of its energy-integrated version one can use the small-cone 
approximation (SGA), as suggested and derived in [1]. The dependence of the impact 
factor on the jet radius R has the form [S] (/x is a shorthand for p/j, p/j and 1 /^) 

dlnu [xj, log R^j = AAxj) \og{R) + (^xj, log ^) + (23) 

and the analytic expressions for the coefficients A, B were explicitly computed |1] for a 
particular jet algorithm (FIP). 

However, whereas the coefficient A of log(i?) depends only on the incoming hadron 
and is given in terms of the usual splitting functions as 

= + Y. At) 

+ [P„(0+2n^P„(C)]/,(^)}. (25) 
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the constant term B depends also on the details of the jet algorithm. Ivanov and Papa 0 
computed such coefficient for an algorithm which was used in pioneering work on QCD 
jets by Furman [7] — we shall refer to it as FIP algorithm — which, however, is no more 
used in present day phenomenology. 

The main purpose of our paper is to derive such coefficient for the two mostly used 
algorithms of QCD analysis, namely the cone-algorithm and the kt-algorithm. The com¬ 
putations can be repeated by following the procedure of |1] . 

The expressions for the jet vertices at LL and NLL level are extracted from the per¬ 
turbative calculation of processes with two incoming partons producing 2-jet at LO and 
NLO respectively. At LO the amplitudes have just two partons in the hnal states, each of 
which is identihed with a jet. The jets are emitted back-to-back in the azimuthal direction 
and have no substructure. Therefore no dependence on the algorithm is found at LO. 

The same is true at NLO as far as the virtual corrections are concerned. On the other 
hand, the NLO real corrections involve 3 partons in the hnal state, therefore a jet can 
be constituted by either one or two of them. In the case of 1-parton (simple) jet, all 
algorithms are designed in such a way that no further emission is found within a region 
of radius R in the (y, 0) plane around the position of that parton. 

Therefore, the differences among the algorithms are to be found in the 2-parton (com¬ 
posite) jet conhgurations. In the following section we will carefully compare the dehni- 
tions the cone, kt, and FIP jet algorithms and, from their differences, we shall compute 
the small-cone impact factors for the cone and kt jets. 


3 Jet algorithms 

3.1 The cone algorithm 

According to ref. [6], when two partons pi and p 2 are combined into one jet of radius R, 
the resulting jet variables {y, E, 0) are dehned to be 


E = El+ E 2 , y = 


yiEi + y2E2 


(j)lEi -f 02^/2 

E 


(26) 


To determine whether the two partons are to be combined, we see if they £t in a cone 
of radius R about the jet axis in the {y, 0) plane. In practice, by denoting with Qij the 
{y, 0)-distance 

^ij = (Vi - Vjf + (0* - (t)jf (27) 

one requires Du < R and D 2 J < R for the composite jet, which amounts the condition 

El -|- E 2 


Di2 < R- 


(28) 


^max(i 7 i, 7^2) 

A simple jet can then be dehned only if the two partons cannot be combined, i.e., provided 

El -|- E 2 


Dio > R 


max(i7i,i72) 


(29) 






3.2 The kt algorithm 

According to ref. [5], the kt clustering algoritm consists in an iterative procedure which 
is based on comparing a set of resolution variables of single-particle dis (B for beam) and 
pairs dij 

diB = E‘f , dij = min(£;2, £;|) ^ . (30) 

One then considers the smallest one: if it is a diB, particle i is thrown in the beam 
basket and removed from the list; if it is a dij then particles i and j are merged into a 
pseudoparticle {ij} — e.g., by using the recombination scheme ([26]). The procedure is 
repeated from the beginning, until all resolution variables are greater than some (hard) 
stopping parameter dent ^ ^qcd- 

In the case of three partons in the hnal state there are these possibilities: 

• One of the diB, say diB is the smallest resolution variable. 

— If diB < dcut then particle 1 belongs to the beam and we are left with two 
partons which, if not in the beam, form two simple jets; 

— If diB > dent then the clustering stops and all three particles form simple jets; 

• One of the dij, say di 2 is the smallest resolution variable. In this case 1 and 2 
are merged into a pseudoparticle {12} and one considers three resolution variables: 
d{i2}B, dsB and d{i2}3- 

— If the smallest is larger than dcut clustering stops and we have one simple jet 
{3} and one composite jet {12}; 

— If the smallest is less than dcut, a (pseudo)particle belongs to the beam and we 
cannot have two jets in the hnal state; 

To summarize, we have a composite jet, say {12}, only if 

di2 < diB Vz Oi2 < R . (31) 

In the other case 

Oi2 > R (32) 

only simple jets are present. 

3.3 The Furman algorithm 

In [1] Ivanov and Papa dehne a jet as a set of particles within a cone of radius R. To be 
more precise, they require that: 

• the jet’s momentum is the sum of the particles’ momenta; 

• all and only the particles of the jet belong to a circle of radius R in the {y, 0) plane 
and centered at the jet’s momentum. 
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This is just the jet definition of Furman [7] used in the ’80s for early phenomenology of 
NLO QCD jets — we denote it “FIP algorithm”. 

In the case under study, where at most two particles (1 and 2) can form a jet, the 
condition for a composite jet is nothing but the prescription adopted in the cone algorithm: 
two particles such that 




< R 


El + E2 
max{Ei, E 2 ) 


(33) 


can be considered as a composite jet. On the other hand, particle 1 can form a simple jet 
if no other particle is found within a distance R from it, namely 


r2i2 > R ■ 


(34) 


This definition is somewhat pathological, because it may happen that a given configura¬ 
tion can give rise to both a composite jet and two simple jets. This fact can be easily 
understood in the case of two particles with the same transverse energy and whose dis¬ 
tance satisfy R < ^12 < 2R. They can form simple jets because a cone of radius R 
centered on either particle does not contain the other; on the other hand, a cone of radius 
R centered halfway the two particles contains both of them. For this reason, it is not 
possible to extend such jet definition into an IR-safe algorithm to all orders, i.e., with an 
arbitrary number of particles, hence it has been abandoned in favour of better algorithms 
like the cone and especially the kt one. 

In any respect, this definition is different from both the cone algorithm (in the case 
of simple jets) and the kt algorithm (in the case of composite jets), and yields different 
results when used to defined any observable. It is the main purpose of this paper to derive 
the correct expressions for the jet vertices in the SCA for the kt algorithm and also for 
the cone one. 


4 Small cone jet vertices 

In this section we shall compute the SCA jet vertex in the cone and kt algoritm. The idea 
is to identify and calculate the contributions differing from those of the original paper [3] 
where the FIP algorithm was adopted. 

As explained in the previous section, for dijet production at NLO, the differences 
among jet algorithms occur only in the way that two partons can be combined to form 
a jet. In the BFKL approach, two partons (1 and 2) in the fragmentation region of an 
incoming hadron are produced by the interaction of a parton stemming from the hadron 
and a reggeized gluon, as depicted in fig. [2j 

In the high-energy kinematics, the incoming parton has just a fraction x of the pure 
longitudinal momentum of the hadron, while the reggeized gluon’s momentum q is essen¬ 
tially transverse. Following Ivanov-Papa (IP) |3] we work in dimensional regularization 
[D = 4 + 2e) and we indicate with ki (/C 2 ) the (2 -|- 2e)-dimensional transverse momentum 
of the outgoing parton 1 (2), and with ( {( = 1 — Q its longitudinal momentum fraction 
with respect to the incoming parton. In terms of these variables, the relative rapidity and 
azimuthal angle between the two outgoing partons are 

= - log , A0 = arccos , C = 1 - C ■ (35) 

2 \ki\ fc2 
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Figure 2: Kinematics of the fragmentation region of hadron A: the collision of the incoming 
parton (blue) and of the Regge gluon (red) produces a pair of outgoing partons (black). In 
parentheses the longitudinal momentum fraction and the transverse momentum of each 
particle. 


We shall present our results by using the notations of ref. |1] , namely in term of their 
quantities, which are related to our dehnition of jet impact factor flTB]) by 


.* / _/V^ — 1 X T 

^nu{y, E, 0 ) = ^ ^ 


1 = Y.K, + h/ 


(36) 


*./ 


where i (/) are labels for the initial (final) state of the sub-processes contributing to the 
cross section, while the superscripts R and V denote real and virtual parts respectively. 


4.1 Vertex for the cone algorithm 

The condition for composite jet in the cone algorithm fl28l) coincides with the one adopted 
by FIP fl5^ . thus no difference is expected for the corresponding contributions. 

In contrast, the conditions for simple jets and fl5T)) are different, causing different 
contributions to the jet vertices. By using the notations of ref. [1], fc = fci is the transverse 
momentum of the parton forming the simple jet, q — k = k 2 is the transverse momentum 
of the parton outside the jet (also called “spectator”), and one introduces the transverse 
vector A 

g = ^ + A (37) 

which vanishes when the two partons are collinear. In fact, for small A, we have 

= Mu + Mu - 

and also 

El + E 2 _ I fell + 1 ^ 2 ! ^ C) _ 1 /ggX 

max(Ei,E 2 ) max(|fci|,|fc 2 |) “ |)t|max(^l,0 max(C, C) 

Therefore, the simple-jet condition 0281) becomes 


C max(C,C) 


(40) 


at variance with the FIP condition fl3^ which, expressed in terms of A, reads like eq. (HOD 
but without the denominator max((C,(^) [Ij. 
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In order to identify the contributions of the FIP jet vertex that have to be modihed 
in the cone algorithm, let us recall that the simple jet conhgurations were computed in 
two steps: by allowing the spectator parton to span the whole phase space and then 2 ) 

by subtracting the contribution stemming from the spectator parton inside the jet cone. 
Therefore, we have to modify only the subtractions by replacing R —)■ R/ max(C, 

The quark initiated jet vertex (quark+Regge-gluon —)■ quark+gluon, see £g. [2]) has 
two such subtractions, one for the quark-jet and one for the gluon-jet, whose results 
are reported in eqs. (5.36) and (5.38) of ref. [1] respectively. In the gluon-jet term the 
substitution R ^ R/ max((C, () is straightforward: 


= (5.36) [1] ^ 


X 


Q, r(i-6) r^(i + e) 
27r e(47r)® r(l -|- 2e) 

R ^ " 


{k 


2N7+6-; 


(k ■ ir f ^ 5 ^ /, 

^ a=q,q 


C 


max(C,C) 


c 


PgqiO ( 1 + 2e log - ) -|- eCpC 


a. 


271^ ’ 


{k-iy 


[ ^ c /a ^OT(C)21og(max(C,C)) , 

^ a=q,q V 8 / 


(41) 


where 7 = ii/ — ^ and I = ei + 162 is a complex vector lying only in the first two of the 
2 -|- 2e transverse dimensions. 

In the quark-jet one has to proceed more carefully, because of the presence of a double 
pole in e multiplying [R/ max((C, ^))^. This would generate, among other things, modihed 
simple poles and also hnite double logs. However, since the double pole is multiplied by a 
5(1 — C), for these terms max((^, C) = 1 ^ind the outcome is identical to the FIP algorithm. 
The only difference comes from the simple pole in front of the Pqq splitting function, and 
we obtain 


= (5.38) [1] 

^ f T ^ (t) ■ 

’ a=q,q V S / 

(42) 

The gluon initiated jet vertex has two subtractions too, one for the qq final state 
and one for the gg hnal state, whose results are reported in eqs. (5.49)0 and (5.56) of 
ref. |1] respectively. The situation is very similar to that of the quark initiated vertex: 
the subtraction in the (anti)quark jet contains a simple pole times the Pqg{C) splitting 
function in front of the R^^ factor, while the gluon jet contains both single and double 
poles, the former with the Pgg{C) splitting function and the latter with the 5(1 — C) 
distribution. For both types of jets the substitution R ^ R/ max((C, () simply amounts 

"‘We note a misprint in ref. [1]: in eqs. (5.48-49) the first subscript of should be g instead of q. Also 
the last subscript in eq. (5.38) should be —g instead of —q. 
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to finite contributions proportional to log (niax((', : 


+ + lf,g-g - (5-49) [l] + {q ^ q} + (5.56) [l] 


jn jK 

g;q,—q ' g-,q,—q ' 9,9,—9 


tR 


tR 


( 43 ) 






(k ■ IT t ? 

a=q,q 




Xj 


[2n/Pgg(C) + Pgg(C)] 2 logniax(C, C) 


To sum up, the jet vertex for the cone algorithm in the small-cone approximation is 
obtained by replacing R ^ R/ max((C, ^) in the hnal formulae (5.39) and (5.57) of |1]. 
The complete expressions are written in eqs. 0531154^ . 


4.2 Vertex for the kt algorithm 

The condition for simple jet in the kt algorithm fl3^ coincides with the one adopted 
by FIP fl3T)l . while the conditions for composite jet fl3T|l and fl3^ are different. In the 
composite jet conhguration the jet’s transverse momentum is k = ki + k 2 = q and it is 
convenient to dehne the auxiliary transverse vector A as 


ki = (k + A 

thus obtaining 

12 ^^ 2^2 

The composite jet condition fl3lD becomes 

|A| <CC|fc|^, 


(44) 

(45) 

(46) 


at variance with the FIP condition fl5^ which reads [U |A| < min(C,C)|fc|^- 

The corresponding contributions of [1], that have to be modihed in order to recover 
the jet vertex in the kt algorithm, are found in sec. 5.1.1.C for the quark initiated vertex 
and secs. 5.2.l.b and 5.2.2.b for the gluon initiated one. In all such cases the modihcation 
amounts to replace 

|Amax| = min(C,C)|fc|^-t CCI^I^ , (47) 

and hnally to substitute min((C, () (( in the (^-integral with the relevant splitting 

function. 

Explicitly, in the quark-initiated case, such integrals for the FIP and kt algoritms 
reads respectively 


jR (FIP) 
^9,9+9 


OC 


jR (kt) 
^9',9+9 


OC 



dC [min(C, C)] 

dC (CC)^^ 


2e l + C^ + <^ 
C 

1 + C^ + 6C^ 
C 



(48) 

(49) 


Their difference (3 — 7r^/3 — 31og2)e, when multiplied by the overall 1/e pole, provides 
a hnite contribution that has to be added to the IP result in order to obtain the proper 
expression for the kt algorithm: 


4= (5.33)[1] ^ ■')" E (3 - y - 31og2) , (50) 

a=9,q ^ 2 
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In the gluon-initiated case the procedure is identical; here we have two contributions: 
one from the qq jet 

4+,- = (5-47)[l] ^ + ^{ky-'Hk . 2>vT„ (I log2 - , (61) 

and one from the gg jet 

= (5-54)|ll ^ ^ “ log2) . 

(52) 


The integrals (1481149^ and the analogous ones for the gluon-initiated contributions, yielding 
eqs. fl5nil52p . were already considered and computed [ 8 ] in the hrst study of the relation 
between the kt and the cone algorithm — the latter sharing with FIP the condition of 
composite jet. 


4.3 Final expressions of the jet vertices 

The result of the jet vertex for the cone and kt algoritms is reported below, by adding to 
the original expressions of IP [1] the modihcations computed in the previous subsections 
and here highlighted in boldface: (• • • )c for the cone and ■ •)k for the kt. The quark 
part is 




P„{C) + ^P „(0 

Uf 


fog—+ 

k-F 


- 2 C-^^[p,,(C) + P,,(C)]fog 


R 


^01 An 

“ Trfog^<^(l - C) 


(max(C,C)>c 2 /ifj 


So 85 TT 


+ {i+e)ic 


P 18 ' 2 ' 2 L 
2 ( 1-0 + 


+ - 0 <1 xl'i] fog ^ + TF + ^ + F V'' (1 + 7 + ?) - 0 (TT - 7 ) - ^ 


n 






C',-& 


C r 21ogC , ^ fiog{l - 0 


+ 5 ( 1-0 


TT^ 9 


TT 


C, ( 31og2 - _ - 1 + ^3 - - - 3 log2)^ j - 


+ CaC + CfC + 


1 + 0 


CaR + 2Ca log I + - 2 logo 


(53) 


where 7 = iz/ — 1 / 2 , (3q = (IIC^ — 4:nfTR)/3, ip is the digamma function, while the 
splitting functions PabiC) and the special functions Ij{n, 7 , () are reported in app. 0 
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The gluon part is 


-2C-""lP„(C) + 2n/P„(C)|log 


P»(C) + ^2nfF„(C) 

Up 


R 


log—+ 

h-F 


_ _ AN 

(max(C, 0 >c ^ 

+ Ca^(1 - C) <! Xnu log ^ ^ k' + 7 + - 7^ - 


£o 1 

P ^ 2 

1 

— T — 

12 6 


Xn 


131 TT^ 11 

-log 2 

36 3 3 /K 


+ 2Ca{1-C^^) 


+ C*A 


i-2 + cc)iogC + ^^Y^^ 


2 + CC 


(i + r'^)x®-2iogc + ^/2 


LC (1-0- 

2n/TR 2^CC + (C^ + C^) + ^-^3^ 

/ 1 /2 23 

+^(l-C)(-Y5 + (3log2-- ^ 


(54) 


It is apparent from the above expression that the log(-R) coefficient A of the jet vertex 
is independent of the jet algorithm, while the constant coefficient B depends on it. 


5 Numerical study 

In this section we assess the quantitative difference among the jet vertices in the three 
algorithms (cone, kt, FIP) that we considered, and also the corresponding accuracy of their 
small-cone approximations (SCA). We shall use the term exact in the sense of “without 
SCA”. 

5.1 Jet vertices versus R 

We start by evaluating the “exact” jet vertices in the three algorithms — we employ the 
numerical code used in ref. [15] — for various values of the jet radius R, and we compare 
them with their SCA. We expect the SCA to be better, the smaller the values of R, 
and increasing discrepancy with increasing R. On the other hand, the differences among 
different algorithms shouldn’t vanish with R, according to our analysis. 

This is actually the case, as can be seen in fig. (S] where we plot the exact (points) and 
small-cone (lines) NLO part of the jet vertex versus R, in the three algorithms mentioned 
before (for a given choice of the parameters n, u, E, y). The common slope of the lines 
represents the coefficient A of log(i?), which is the same for the three algorithms, while 
the intercepts at i? = 1 give the constant coefficients B, which clearly depend on the 
algorithm. A detailed study, carried out with several values of the parameters, shows 
that the small-cone approximation works very well up to R of few tenths, with an error 
below 1% for R < 0.2 which increases up to 3-4% when R = 0.5. 
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Figure 3: Dependence on the jet radius R of the NLO jet vertices in the three jet 
algorithms discussed in the text: FIP (dashed red), cone (dotted green) and kt (solid blue). 
The big dots correspond to the exact evaluation (the small error bars showing MonteCarlo 
integration uncertainties), while the straight lines denote the small-cone approximation. 
Here n = D, v = D, y = 3.6, E = 35 GeV. 

5.2 SCA versus algorithm choice 

Next, we specialize our analysis to realistic values of the jet radius and energy. Since 
the typical phenomenological studies on MN jets use i? ~ 0.5 and jet transverse energies 
E >35 GeV, in the following all quantities will be evaluated at R = 0.5 and E = 35 GeV. 
In addition, nowadays the mostly used jet algorithm is the kt. Therefore we adopt “exact” 
quantities, computed in the kt algorithm, as reference quantities, and we estimate the 
deviations to them introduced by the SGA. 

To some extent, adopting the SGA at fixed i? is a sort of choosing a jet algorithm. 
The natural question then arises: how does the discrepancy introduced by the small-cone 
approximation 

(exact-kt) — (SGA-kt) (55) 

compare with the discrepancy caused by different choices of jet algorithm, i.e., 

(SGA-kt) - (SGA-FIP) ? (56) 

In order to answer this question, we compute the exact NLO part@ of jet vertex 
in the kt algorithm as function of z/, and compare it with the SGA in the same algorithm 
and also with the SGA in the FIP algorithm. 

Fig. m shows such a comparison for n = 1. We can see that the exact result in the kt 
algorithm is much better approximated by the SGA in the same algorithm rather than by 
the FIP choice, in particular at small values of u, which are the most important in the u 
integrals of eqs. fl20ll22jl — since the GGF is peaked around z/ = 0. 

This conclusion is further supported by analysing the whole (LO-I-NLO) z/-integrand of 
eq. fl2^ . again by comparing exact-kt, SGA-kt and SGA-FIP, as in hg. |5l The discrepancy 

®We recall that the LO part is independent of the jet algorithm. 
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Figure 4: Comparison of the NLO exact jet vertex in the kt algorithm (points) with its 
SC A in the same algorithm (solid blue) and in the FIP algorithm (dashed red). Here 
R = 0.5, n = 1, y = 3.6, E = 35 GeV. 

introduced by the SCA in the wrong FIP algorithm is about three times larger than that 
introduced by the SCA in the proper kt algorithm, the latter being of the order of 5%. 

Actually, the relative error due to the SCA is slightly larger for the full (LO+NLO) 
quantities than for the pure NLO ones. This is due to the fact that the NLO corrections 
usually have sign opposite to the common LO terms, giving rise to cancellations in their 
sum that amplify the relative differences. 

5.3 Cross section and angular coefficients 

Finally, we present the results of the differential cross section da/AY and of few angular 
coefficients Cm/Cn- 

In £g. [H]we plot the differential cross section da /AY = Cq{Y) by comparing again the 
exact kt calculation with the small-cone approximations in the KT and FIP algorithms. 

It is evident that the wrong choice of the algorithm yields a large error, especially at 
lower values of K, while the sole SCA with the proper algorithm introduces an error of 
4-8%. 

The shape of the Cq curves, which are not monothonically decreasing in Y as one 
could naively expect, is due to an additional cut in rapidity \yi\ > jV min = 3 that we have 
imposed just for computational convenience, as will be shortly explained. Due to this 
cut, the minimum value of Y that we allow is Kmin = = 6 and in this limit the 

cross-section vanishes. It then quickly rises for Y > 6 before eventually decreasing at 
larger Y >7. 

The reason for imposing the ly) > ymin cut is due to the fact that monochromatic (hxed 
E) impact factors oscillate at large z/, causing the u integration in eq. fl^ to be slowly 
convergent. Actually, the convergence is provided by the GGF, whose modulus decreases 
at large u, the decrease being faster at larger values of Y. On the other hand, at low values 
of Y, and in particular with very asymmetric rapidity conhgurations (|||/i| — \y 2 \\ Y 1), 
the numeric z/-integration has to be pushed to large z/-values, thus demanding a large 
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Figure 5: Comparison of the exact v-integrand for kt algorithm (points) with the SC A 
in the same algorithm (solid blue) and in the FIP algorithm (dashed red). On the left: 
n = 0; on the right: n = 1. The upper plots display absolute values in arbitrary units, 
and below them we show the ratios w.r.t. exact integrand. The parameters are R = 0.5, 
El = E 2 = 35 GeV, yi = 3.6, j /2 = 2.8. 


computational effort. Since this problem is absent for realistic phenomenological studies 
(where the impact factors are integrated in the E variable and decrease themselves with 
z/), and because the goal of this analysis is just to compare the main features of the 
different algorithms, we solve the convergence issue by imposing the mentioned cut in 
rapidity \yi\ >3. 

From the plots of figs. IHII^ one can infer that the main effect of the SCA is mostly 
an overall normalization change, and that this is also true even for the wrong choice of 
algorithm (kt versus FIP), though with a larger factor. If this were the case, by comput¬ 
ing ratios of observables such effects should cancel out and reproduce more faithfully the 
exact quantities. This is partially true, as can be seen in £g. [TJ where we plot some ra¬ 
tios of angular coefficients Gm(T)/G„(F) = (cos(mA0j))/(cos(uA0j)) which are usually 
adopted in order to measure the azimuthal decorrelation of the MN jets. 

It is nevertheless evident that the cancellation of the systematic effects is more effective 
if the SCA is made with the proper algorithm, leading to a discrepancy of about 2% or 
less for all the ratios considered. A different choice of algorithm yields definitely larger 
discrepancies, and therefore should be avoided. 


6 Conclusions 

In this paper we have reconsidered the high-energy factorization formula for Mueller- 
Navelet jet production, and we have computed the NLO jet vertices in the small-cone 
approximation for two different jet algorithms: the kt algorithm and the cone algorithm. 
The small-cone approximation amounts to evaluate the generic jet vertex T at small 
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Figure 6: Comparison of the exact differential cross section da/dY (n = d) for kt algo¬ 
rithm (solid blue) with the SCA in the same algorithm (dash-dotted blue) and in the FIP 
algorithm (dashed red). Top: absolute values in linear scale; bottom: ratios of the SCAs 
w.r.t. the exact one. Here R = 0.5, Ei = E 2 = 35 GeV. 

jet radius R ^ 0 according to the expansion \I/ = y41og(i?) + i? + 0{R^), where the 
coefficient A is universal, depending only on the infra-red properties of QCD, while the 
constant term B depends on the jet algorithm. It should be noted that the neglected 
term in the expansion is quadratic in R, and that the linear term is missing. It turns out 
that the small-cone approximation is quite accurate for realistic values of 7? ;< 0.5 and, 
because of its analytic and computational simplicity, it has been and will be used with 
beneht. 

The calculation of the Mellin-Fourier components of the jet vertices (cfr. sec. 12.2p 
was originally performed |1] by adopting an algorithm which is equivalent to Furman’s 
one [7], but which is not used anymore in modern phenomenological studies, being infra¬ 
red unsafe. In order to use such vertices for the analysis of data from high-energy colliders, 
we have computed them in the two most popular algorithms: the kt and the cone ones. 
The ensuing analytical expressions have been taken out in sec. 01 and the differences 
among the three algorithms have been highlighted. 

The quantitative difference between the small-cone approximations in the kt algorithm 
and in Furman’s one have been estimated in sec. |5], by plotting several quantities of 
interest: the jet vertices, the integrand of the high-energy factorization formula, the 
differential cross section w.r.t. the rapidity distance Y = ||/i — 1 / 2 I of the two MN jets, and 
some angular coefficients indicating the azimuthal decorrelations of the jets themselves. 
Also the corresponding quantities with the exact jet vertices in the kt algorithm have 
been plotted. 
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Figure 7: Comparison of some ratios Cm/Cn determined with the exact jet vertices in the 
kt algorithm (solid blue) with those obtained in the SCA with the same algorithm (dotted 
blue) and with the FIP algorithm (dashed red). Here R = 0.5, Ei = E 2 = 35 GeV. 

It turns out that the difference between the kt and Furman’s algorithms is sizeable, 
of the order of 20% at the level of cross section, and about 5% for the ratios Cm/Co = 
{cos{mA(()j)) and Cm/Cn, the angular coefficients of azimuthal decorrelation. 

On the other hand, the discrepancies between the exact results and the small-cone 
approximated ones are much smaller, of the order of 5% at the level of cross section, 
and less than 2% for the angular ratios. We therefore conclude that the small-cone 
expansion, computed with the proper jet algorithm, provides a good approximation to 
the Mueller-Navelet jet vertices; it can then be used as a very convenient tool to perform 
phenomenological studies, in that it requires much less computational resources than the 
exact computation. 

This aspect could be essential when analysing observables obtained by integrating the 
jet energies in a non-factorized domain, e.g., by requiring E 1 + E 2 > 2Ecnt- Such condition 
is often used in dijet analysis since, at fixed NLO, it yields more stable results than the 
condition Ei > E mm , while retaining the 1 y-)- 2 symmetry. On the other hand, a non- 
factorizable domain of integration in energy prevents the use of eq. fl22l) with independently 
integrated impact factors ; a numerical integration in energy (and possibly in rapidity) 
of the integrand fl2UD is thus necessary, and the small-cone approximation represents a 
valuable tool to reduce the computational effort of such calculation. 
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A Expressions of splitting and special functions 

The splitting fnnctions fonnd in the jet vertices are defined in the usnal way: 

'I + z^' 


Pqq{^) — 


Pgqi.^) — Cp 


l-z ^ , 

1 + (1 - z} 


= c. 


I + I 3 ^ ^ 

L(i-^)+ 2 


P,,{z) = Tn [z^ + (1 - 

Pgg{z) = 2 Ca 


11 . ^ 
H-2 + z[l — z) 


_(l-2;)+ ^ 


(57a) 

(57b) 
(57c) 

+ I ^0.4 - ^) i(l - . (57d) 


The fnnctions Ij are expressed in terms of hypergeometric fnnctions and read [ 1 ] 


^2(n,7,C) = I 


c 


2^1(1 ,1 + 7 — 2 + 7 — ^, ((■) 27^1(1 ,1 + 7 + 2 + 7 + |, ((■) 


t-7-1 f+7+1 

^ 2 ^ 1 ( 1 , -7 - f, 1 - 7 - f, 0 2^1(1, -7 + t, 1 - 7 + f , 0 


+ C 
+ (1 + C-^^) 


7i,3(’^, 7, C) = 7, C) ± ^ ■ 


t + 7 

■21ogC) +21ogC 
1 - 


t-7 


logC + 


and we recall the definition 7 = iz/ — 1/2. Xnu is dehned in eq. (ITT)) . 
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(0) 


- 21 og^ 
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